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! Based on cosmological rates, it is probable that at least once in the last Gy 

the Earth has been irradiated by a gamma-ray burst in our Galaxy from within 
2 kpc. Using a two-dimensional atmospheric model we have performed the first 
computation of the effects upon the Earth's atmosphere of one such impulsive 
! event. A ten second burst delivering 100 kJ/m 2 to the Earth penetrates to 

the stratosphere and results in globally averaged ozone depletion of 35%, with 
depletion reaching 55% at some latitudes. Significant global depletion persists 
for over 5 years after the burst. This depletion would have dramatic implications 
for life since a 50% decrease in ozone column density results in approximately 
^ \ three times the normal UVB flux. Widespread extinctions are likely, based on 

extrapolation from UVB sensitivity of modern organisms. Additional effects 
include a shot of nitrate fertilizer and NO2 opacity in the visible providing a 
cooling perturbation to the climate over a similar timescale. These results lend 
^ I support to the hypothesis that a GRB may have initiated the late Ordovician 

mass extinction (Melott et al. 2004). 
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1. Introduction 

Gamma-ray bursts (GRB) within our galaxy have been suggested as a possible threat 
to life on Earth (Thorsett 1995; Scalo & Wheeler 2002; Dar & De Rujula 2002; Melott et al. 
2004). Some effects similar to that due to a nearby supernova (SN) are expected, especially 
depletion of stratospheric ozone due to ionization caused by incident gamma radiation. GRB 
are rarer than supernovae, but their greater energy output results in a larger region of 
influence and hence they may pose a greater threat. It is likely (Melott et al. 2004; Dermer 
& Holmes 2005) that in the last Gy a GRB has occurred close enough to have dramatic 
effects on stratospheric ozone, leading to detrimental effects on life through increases in 
solar ultraviolet (UV) radiation which is strongly absorbed by ozone. A major issue is the 
timescale for atmospheric chemistry: most of the GRB fluence comes in seconds or minutes 
versus months for supernovae. 

In order to gain more detailed and accurate insight into these expected effects we have 
performed computations using the Goddard Space Flight Center (GSFC) two-dimensional 
atmospheric model. This model has been used previously to investigate the atmospheric 
effects of SN (Gehrels et al. 2003). While we build on that work, the computations discussed 
here are significantly more challenging due to the extremely short duration and greater 
energy output of GRB in comparison to SN. 

2. Methods 

We take as "typical" a GRB with power 5 x 10 44 W (isotropic equivalent) and duration 10 
s, whose gamma-ray spectrum is described by the Band spectrum (Band et al. 1993). These 
assumptions are drawn from observations and are not dependent upon beaming angle. We 
have determined the depletion of ozone for such a GRB beamed at the Earth from a distance 
of 2 kpc, delivering to the Earth a total fluence of 100 kJ/m 2 . This distance corresponds to 
that of a probable nearest "typical" GRB in the last Gy, based on conservative assumptions 
(Melott et al. 2004). 

The prompt effect of this burst at the Earth's surface is a "flash" of UVB radiation 
with power ~ 20 W/m 2 (Smith et al. 2004). This is about seven times the intensity at the 
Earth's surface on a bright, sunny day, but is brief and so is not likely to have a major effect 
on life. Longer term effects include ozone depletion and the resulting increase in solar UVB 
flux, which we begin to explore here. 

We have not included the effects of any ultra high energy (> 10 18 eV) cosmic rays from 
a GRB (Dermer & Atoyan 2004; Waxman 2004a,b; Dermer & Holmes 2005), due to uncer- 
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tainty as to whether and at what energies GRB may produce such particles. Lower energy 
cosmic rays which may be produced are deflected at our assumed distance by the galactic 
magnetic field. Inclusion of ultra high energy cosmic rays would measurably increase the 
ozone-destroying energy budget. In this case, we would also expect production of some ra- 
dioisotopes from cosmic ray spallation off of atmospheric constituents, as well as a terrestrial 
burst of muons. 



2.1. Atmospheric Model 

The GSFC 2D model is described in Douglass et al. (1989); Jackman et al. (1990); 
and Considine, Douglass, & Jackman (1994). The model's two dimensions are latitude and 
altitude (ranging up to about 116 km). The latitude range is divided into 18 equal bands and 
extends from pole to pole. The altitude range includes 58 evenly spaced logarithmic pressure 
levels (approximately 2 km spacing). A lookup table is used for computation of photolytic 
source term, used in calculations of photodissociation rates of atmospheric constituents by 
sunlight (Jackman et al. 1996). Winds and small scale mixing are included as described in 
Fleming et al. (1999). For this study, we have removed anthropogenic compounds such as 
CFCs. 

We have employed two versions of the atmospheric model. One is intended for long term 
runs (many years) and includes all transport mechanisms (e.g., winds and diffusion); it has 
a time step of one day and computes daily averaged constituent values. The second is used 
for short term runs (a few days) and calculates constituent values throughout the day and 
night, but does not include transport. Previously, this version has been used with a time 
step of 225 seconds (Jackman et al. 2001). In the current study we have used a time step 
of one second in order to allow for spreading our GRB gamma radiation input over several 
time steps. 

Gamma-rays are introduced in the model in a manner similar to that described in 
Gehrels et al. (2003). In that study, gamma-rays were included using the spectrum of 
SN 1987A. In the present study, the gamma-ray differential photon count spectrum used 
is that of Band et al. (1993), which consists of two smoothly connected power laws. We 
use the following typical values for the break energy and power law indices, respectively: 
E = 187.5 keV, a = -0.8, (3 = -2.3 (Preece et al. 2000). The total incident energy is scaled 
to our desired value (in this study, corresponding to a fluence of 100 kJ/m 2 ) The total photon 
flux in each of 66 evenly spaced logarithmic energy bins, ranging 0.001 ^ E ^ 10 MeV, is 
obtained by integrating the Band spectrum for each bin. 
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2.2. Simulations 

All simulation runs used for analysis were begun with initial conditions obtained from 
a long-term (roughly 40 years) run intended to bring the model to equilibrium. Constituent 
values from this run are read in by the Is time step version of the model which runs for 7 days 
(beginning at noon), either with or without input of gamma-rays. Runs including gamma 
radiation treat the burst as a step function at noon on day 4, with duration 10 s. In the 
current study, the burst is input in late March (near the spring equinox) with incident angle 
degrees (equatorial). (Forthcoming studies will investigate the effects of varying intensity, 
incidence angle, and times of year at which the burst occurs.) Day 4 is the middle of the 
short term run, which allows for "warm up" in the model and insures that relevant chemistry 
is accurately computed over many time steps after the burst. Constituent values from this 
type of run are then read in by the 1 day time step version which is run for 20 years in order 
to investigate long term effects and determine how long the atmosphere takes to return to 
equilibrium, pre-burst conditions. 

Analysis is then performed by comparing such a combined base-short-long run without 
gamma-ray input to such a run with the burst included. Ozone depletion and other effects 
are computed by comparing these two runs. 

3. Results 

Stratospheric ozone is lost through several catalytic reactions involving oxygen-, nitrogen- 
, hydrogen-, chlorine-, and bromine-containing gases. The constituents in the stratosphere 
are generally grouped into "families" such as O x (0 3 , O, 0( 1 D)), NO y (N, NO, N0 2 , N0 3 , 
N 2 5 , HN0 3 , H0 2 N0 2 , C10N0 2 , BrON0 2 ), HO x (H, OH, H0 2 ), Cl y (chlorine-containing 
inorganic molecules), and Br y (bromine-containing inorganic molecules), which allow for ef- 
ficient computation of the chemistry and transport effects. We assume that there were no 
anthropogenic sources for any of these families, as human influence has been negligible for 
most of geologic time. 

In the case of a large input of gamma rays to the atmosphere NO y compounds (most 
importantly NO and N0 2 ) are created through dissociation of N 2 in the stratosphere which 
then reacts quickly with 2 to generate NO. Subsequent reactions create N0 2 and other 
compounds. Together, these react catalytically to deplete O3 through the cycle 



NO + 3 — N0 2 + 2 
N0 2 + O -> NO + 2 



(1) 
(2) 
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net : 3 + O -> 2 + 2 (3) 

Note that NO is not consumed in this cycle and the net result is the destruction of O3 
and production of 2 . Other reactions can complicate this cycle, such as destruction of NO 
by reaction with N; production of O3 through reactions of NO with H0 2 ; and interference of 
NO y with other families (chlorine-, bromine-, and hydrogen-containing constituents) which 
reduces the ozone depletion from these families. Some uncertainties in the atmospheric 
model's treatment of this cycle are discussed in Section 4. 

The primary results of our simulations are increases in NO y and decreases in O3. Ozone 
column densities can then be used to calculate the resulting UVB flux at the Earth's surface 
and from this information biological effects can be estimated. UVB is particularly dangerous 
to organisms because DNA is damaged by absorption in this wavelength range. 

Results of our modeling are shown in Figures 1, 2 and 3. We have modeled the effects 
of 100 kJ/m 2 total incident gamma-ray fluence, input as described in Sect. 2.2. This corre- 
sponds to our "typical" GRB located at about 2 kpc. Figure 1 shows the vertical column 
density of NO y at each latitude over time. The burst is input at time 0. Figure 2 shows the 
vertical column density of O3. Included are scales in both Dobson units (the usual unit of 
ozone column density) and 10 18 cm -2 . A Dobson unit describes the thickness of a column of 
ozone at STP and is defined as 1 DU = 0.01 mm thickness (or, 1 DU = 2.69 x 10 18 cm -2 ). 
Figure 3 shows the percent difference at a given location (between a run with gamma-ray 
input and one without) in vertical column density of O3. Immediate depletion of ozone 
is evident. Due to their qualitative similarity, we have chosen not to plot here changes in 
NO y . Maximum increase in NO y is largely coincident with maximum decrease in O3. The 
maximum increase in NO y at a given location is approximately 30-fold for this case. 

Several features in these plots are worth noting. First, as is seen in Fig. 2 and 3, 
depletion of ozone is initially greatest at the equator (where the incident flux is highest), 
becoming greatest toward the poles within a year or so. Larger ozone depletions at the 
poles are primarily due to the long lifetime of the enhanced NO y in the polar stratosphere. 
Figure 3 gives a somewhat exaggerated impression of the effect of depletion at the poles, 
since ozone is initially high there. The enhanced NO y , including HNO3, will lead to an 
enhancement of nitric acid trihydrate (NAT) polar stratospheric clouds (PSCs). These NAT 
PSCs facilitate heterogeneous reactions that result in greater ozone depletion by halogen 
(chlorine and bromine) constituents. This is especially true in the south polar region where 
a stronger polar vortex with colder stratospheric temperatures is in place during the winter. 
This contributes to an asymmetry which would be peculiar to the present-day configuration 
of continents. A much larger effect which contributes to the polar asymmetry is the time of 
year at which the burst occurs, since ozone concentrations at the poles exhibit large seasonal 
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variations. This asymmetry is largely north-south reversed for a burst in September rather 
than March (Thomas et al. 2005, in preparation). We therefore conclude that present 
peculiarities of continent distribution (including the south polar vortex) are not a major 
source of uncertainty in this work. 

Around 5-6 months after the burst there is a short-lived production of ozone toward the 
south pole. Production occurs at the end of south polar night when a lack of photolysis has 
caused accumulation of NO y constituents which are suddenly photolyzed as the sun rises, 
producing O which may then react with O2 to form O3. 

Globally averaged ozone depletion reaches about 35% (at the start of the long-term 
run) and a maximum depletion of about 55% occurs first at the equator immediately after 
the burst, and then again about 15 months after the burst, in the southern hemisphere. 
Significant global depletion (10% or more) lasts for over 5 years after the burst. 

Figure 4 (available as an animation online at http://kusmos.phsx.ukans.edu/~melott 
/DNA_damage_MPEG.mpg) shows DNA damage estimated by convolving the daily average 
UVB flux at the ground with a biological weighting function (Setlow 1974; Smith et al. 1980). 
We include only ozone absorption effects on the UVB flux since the effect of scattering 
at these wavelengths is comparatively small. We have normalized the plot by dividing 
the damage by the annual global average damage in the absence of a GRB. Greater DNA 
damage probability is evident at low latitudes. This is due to combination of the O3 depletion 
effects with the sun incidence angle, length of day, etc. We have performed other runs with 
different GRB incidence latitudes, times of year, etc., to be discussed elsewhere (Thomas 
et al. 2005, in preparation) and found that the concentration of computed DNA damage 
to low-mid latitudes is a general feature of the GRB hypothesis. One might think that 
this damage would be countered by a greater evolved UVB-resistance in organisms at low 
latitudes. However, at least for modern organisms, there is no evidence that temperate zone 
phytoplankton are any more UVB- resistant than Antarctic plankton (Prezelin 2004). Thus, 
one might predict that greater ecological damage and extinction would be likely near the 
equator. It is interesting that the late Ordovician mass extinction seems to be alone in 
having recovering fauna preferentially derived from high-latitude survivors (Jablonski 2004; 
Sheehan 2001). 

4. Uncertainties 

Our ionization profiles are computed using simple energy-dependent attenuation co- 
efficients, instead of a full radiative transfer calculation. This technique is implemented 
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following Gehrels et al. (2003) with the primary modification being the functional form of 
the spectrum. That study found good agreement between their calculations and a full ra- 
diative transfer model. We find that ionization due to the gamma-ray input peaks around 
30 km elevation, which is in agreement with Gehrels et al. (2003) and Smith et al. (2004). 

Effects of prompt redistributed UV are not included in the model. This will both produce 
and destroy ozone. We have performed a simple test of the magnitude of such effects by 
increasing the solar flux by 10 4 times for 10s to simulate the redistributed UV. The resulting 
ozone depletion is a few percent greater. Also, in some cases the ionizing fluence from the 
GRB afterglow may be as large as that from the burst. So, at worst our results are somewhat 
conservative. As discussed in Melott et al. (2004) and Dermer & Holmes (2005), there is 
enough uncertainty in the GRB rate at low redshift that the fluence from the probable 
nearest burst could be one order of magnitude greater than we model here. 

A discussion of some uncertainties in the production of NO y compounds is presented in 
Melott et al. (2004). In particular, reactions involving exited state nitrogen atoms, N( 2 D), are 
not included in the atmospheric model. Comparison with simplified, off-line computations 
using various ratios of N( 2 D)/N( 4 S) indicate that for our present case the effect on NO 
production of including excited state N atoms is small. Though results depend on the 
assumed temperature at which the reactions occur, for a reasonable range of temperatures 
(230-270 K), there is little variation of NO production with increased concentrations of 
excited state N atoms, and there is generally less than a factor of two difference between the 
off-line computations and model results. An additional complication is the effect of reactions 
involving HNO which may limit the production of NO y . Such reactions are not included in 
the model, but we estimate that this is a small effect. 

The GSFC two-dimensional atmospheric model is empirically based and its dynamics 
are not coupled to the significantly changing constituent levels and accompanying heating. 
This fact introduces some uncertainty in the transport of constituents. As mentioned before, 
variations due to input latitude and time of year are likely more significant. 

5. Discussion 

A significant result of our modeling is that even for a short duration input of radiation, 
atmospheric effects are large and long-lived. Expectations based on supernova studies in- 
dicated that a short duration input might not have such large effects (Gehrels et al. 2003). 
Appearance of features such as the localized production of 3 highlights the need for detailed 
modeling of the effects of a GRB on the Earth's atmosphere, as complexities and feedbacks 
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in the chemical processes can be important and simple scaling may not prove reliable. 

Melott et al. (2004) summarizes studies of UVB sensitivity of various organisms (see 
also Cockell (1999)). About 90% of UVB is presently absorbed by atmospheric ozone. Due 
to the sensitivity of DNA to this radiation, increases of only 10-30% can have lethal effects on 
many organisms, especially phytoplankton, the base of the food chain. Ozone depletions in 
the range of 50%, as seen here, lead to roughly three times more UVB at the surface, which 
is clearly a possible candidate for causing mass extinctions. Of course, we expect additional 
events from smaller burst fluences over the last Gy, less intense but still significant for the 
biosphere. 

There are other effects. The event described here could potentially produce of order 
0.5 g/m 2 mean global deposition of nitrates. Biota are generally nitrate-starved, and this 
deposition may have eased the transition to land, which accelerated after the Ordovician. 
This nitrate deposition may provide a geochemical signature which could serve as a test of 
our hypothesis, though this would be difficult due to the extreme water solubility of nitrates. 
On the other hand, our hypothesis is falsifiable on geochemical grounds. That is, a layer of 
iridium (associated with impact events) or radioisotopes such as 244 Pu (associated with SN 
events (Ellis et al. 1996)) would not be associated with our scenario. Different radioisotopes 
could be generated by spallation if significant levels of cosmic rays are received (see Sect. 2). 
However, few would survive to the present for the late Ordovician mass extinction, which 
occurred 443 My ago. 

The Ordovician extinction is associated with a brief glaciation in the middle of a period 
of stable warm climate. We speculate that there may have been a significant perturbation 
by the opacity of NO2, which would cut off a few percent (ranging up to 35% for a month or 
so during polar fall) of solar radiation (Reid et al. 1978). This would occur primarily at high 
latitudes, as can be seen in Fig. 1. The removal of O3 (a greenhouse gas) also may cause 
some cooling, but this effect should be negligible compared to that due to the increase in 
N0 2 We will provide more detail on these ideas in the near future. 
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Fig. 2. — Column density of O3 with scales for both Dobson units (left) and 10 18 cm 
(right). (Plotted from one year before the burst to four years after.) 
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Months 

Fig. 3. — Pointwise percent change in column density of ozone (comparing runs with and 
without burst). 
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Fig. 4. — Relative DNA damage (dimensionless), normalized by the annual global average 
damage in the absence of a GRB. (Plotted from one year before the burst to four years 
after.) Note that white is set to 1.0. This figure is available as an animation online at 
http:/ /kusmos. phsx.ukans.edu/~melott/DNA_damage_MPEG.mpg . 



